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The impact of selection biases on the ^peak — ^iso correlation of 
Gamma Ray Bursts 
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ABSTRACT 

We study the possible effects of selection biases on the -Bpoak — ^iso correlation caused by the 
unavoidable presence of flux-limits in the existing samples of Gamma Ray Bursts (GRBs). 
We consider a well defined complete sample of Swift GRBs and perform Monte Carlo simula- 
tions of the GRB population under different assumptions for their luminosity functions. If we 
assume that there is no correlation between the peak energy i^pcak and the isotropic luminos- 
ity Liso, we are unable to reproduce it as due to the flux limit threshold of the Swift complete 
sample. We can reject the null hypothesis that there is no intrinsic correlation between ii'pcak 
and Liso at more than 2.1a level of confidence. This result is robust against the assumptions of 
our simulations and it is confirmed if we consider, instead of Swift, the trigger threshold of the 
BATSE instrument. Therefore, there must be a physical relation between these two quantities. 
Our simulations seem to exclude, at a lower confidence level of 1.6cr, the possibility that the 
observed ii^pcak-iiso correlation among different bursts is caused by a boundary, i.e. such that 
for any given ii^peak, we see only the largest Liso, which has a flux above the threshold of the 
current instruments. 

Key words: Gamma-ray: bursts — Radiation mechanisms: non thermal 



1 INTRODUCTION 

The isotropic luminosity Liso and the isotropic energy _Biso of 
Gamma Ray Bursts (GRBs) are strongly correlated with their peak 
spectral energy -Bpcak (i-e. the peak of the vF^, spectrum). The 
spectral-energy correlations of GRBs gathered the interest of the 
scientific community for their possible implications on the physics 
of GRBs (Yamazaki, loka & Nakamura 2004; Eichler & Levinson 
2005; Lamb, Donaghy & Graziani 2005; Levinson & Eichler 2005; 
Rees & Meszaros 2005; Toma, Yamazaki & Nakamura 2005; Bar- 
biellini et al. 2006; Ryde et al. 2006; Thompson 2006; Giannios & 
Spruit 2007; Thompson, Meszaros & Rees 2007; Guida, Bemar- 
dini & Bianco 2008; Panaitescu 2009), but also raised an intense 
debate about the strong selection biases they might suffer of (Band 
& Preece 2005; Nakar & Piran 2005; Butler et al. 2007; Butler, 
Kocevski & Bloom 2009; Shahmoradi &Nemiroff 2011; Kocevski 
2012; but see Ghirlanda, Ghisellini & Firmani 2005; Bosnjak et al. 
2008; Ghirlanda et al. 2008; Nava et al. 2008; Amati, Frontera & 
Guidorzi 2009; Krimm et al. 2009). 

Ghirlanda et al. (2012) studied the comoving frame properties 
of a small sample of GRBs for which the bulk Lorentz factor F 
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could be estimated. They found that F is correlated with L'f^^ and 
with -Epcak- These newly found correlations offer an interpretation 
of the -Epoak — Liso coiTelation as a sequence of F factors, the most 
luminous GRBs have the largest F and Lpeak- 

Since these correlations involve the rest frame prompt 7-ray 
emission properties (Lpcak, Liso and Liso) of GRBs with measured 
redshift z, it has been argued that they are strongly biased by (i) the 
instrumental detector-trigger threshold and/or by (ii) the redshift 
dependence of the involved observables (e.g. Butler et al. 2007). 

Different studies (Ghirlanda et al. 2008; Nava et al. 2008) 
quantified the possible instrumental selection biases finding that, 
even if present, they cannot be responsible for the existence of the 
spectral-energy correlations. While the Lpcak — Liso and Lpoak ^ 
Liso correlations are defined considering the time-integrated spec- 
tra of GRBs, it has been shown (Firmani et al. 2009; Ghirlanda et 
al. 2010; 2011a; 2011b) that similar correlations hold within indi- 
vidual bursts. This is a strong argument in favour of the physical 
origin of the correlations since within a single burst the instrumen- 
tal selection effects (e.g. the flux-limited threshold of any burst de- 
tector), or possible effects due to the dependence of Lpcak, Liso 
and Liso on the redshift do not play a role. 

Discovered with a dozen of GRBs with known redshifts (Yo- 
netoku et al. 2004; Amati et al. 2002), the Lpoak - Liso and the 
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Epeak — Eiso coiTelations have been updated in the last decade by 
an increasing number of bursts with measured z (see e.g. Nava et al. 
201 1 for a recent update). Nevertheless, these samples are not com- 
plete in the sense that there is no well defined flux limited sample 
with the redshifts measured for all GRBs. This raised the suspect 
that the incompleteness in the redshift knowledge might bias such 
correlations. 

Nava et al. 2012 (N12, hereafter) studied the Speak - -Eiso 
and the Bpeak — ^iso correlations with a flux limited sample of 
GRBs detected by Swift which has also a high level of redshift de- 
termination. This complete Swift sample, presented in Salvaterra 
et al. 2012 (S12, hereafter), was selected starting from the 50% 
redshift complete sample of Jackobssorjj and considering the 58 
bright Swift GRBs with a 1 s peak flux (integrated in the 15-150 
keV energy band) F ^ Fum = 2.6 ph cm^^ s~^. The redshift 
recovery rate of a such selected sample turns out to be 90%. There- 
fore, this flux-limited sample has also a high level of complete- 
ness in redshift. S12 study the luminosity function (LF) of GRBs 
throughout the Swift complete sample finding that evolution in lu- 
minosity or in density is required in order to account for the obser- 
vations. Among other results, N12 find that the correlations defined 
with the complete Swift sample are statistically robust: the rank 
correlation coefficient is p = 0.76(0.70) and chance probability 
P = 7 X 10-1° (1 X 10-«) for the Sp^ak - £^1.0 (Sp^ak - iiso) 
correlation. Moreover, the correlation properties (slope and normal- 
ization) of the complete Swift sample are consistent with those de- 
fined with the incomplete larger sample of 136 bursts with known 
z and spectral parameters (see N12). 

Complete samples are at the base of any population studies. 
In this paper we study the impact of instrumental selection effects 
on the -Epcak ^ ^iso coiTelation using as a reference the complete 
Swift sample of GRBs described in S12 and analyzed (in terms of 
its -Epcak — Lisa Correlation) in N12. In particular, since no clear 
consensus has been reached yet on the physical origin of these 
correlations we aim to answer to this specific question: might the 
Epcsk — Lisa correlation be produced by the threshold of a flux- 
limited sample of bursts? 

To this aim we use a population synthesis code that generates 
a large sample of GRBs (following some prescriptions for its lu- 
minosity and redshift distribution - described in §2) and assuming 
that there is no correlation between -Epcak and Lisa- Throughout 
the paper we make use of the following nomenclature: 

• "complete Swift sample": this is the sample defined in S 12 of 
58 GRBs detected by Swift-BAT with F > Fun„ of which 90% 
have measured redshifts. Forty-six out of 58 bursts (79%) have 
well constrained peak energy and Liso and their -Epcak — ^iso cor- 
relation (which is the term of comparison of our simulations) has 
been published in N12. The other 12 GRBs are consistent with the 
£^pcak — Liso correlation as discussed in N12. 

• "complete sample of simulated bursts": this is the sub-sample 
of simulated GRBs with F ^Fum, i.e. representative of the "com- 
plete Swift sample". 

By applying the Swift flux limit Fum (i.e. the same used to define 
the complete Swift sample in S12), we test if a strong Speak ^ 
Liso correlation is obtained in the complete sample of simulated 
bursts (§3) as only due to the cut in the flux. If we cannot find 
a statistically significant -Epeak ^ Liso correlation in the simulated 
complete sample we then reject the null hypothesis that the real one 
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is due to the flux-limit of the sample selection. In the latter case, we 
are left with two possibilities, that we also test in this work (§5), for 
the nature of the Epeak — Liso correlation: either it is an intrinsic 
correlation to GRBs (as also supported by the existence of a similar 
correlation within individual bursts - Ghirlanda et al. 2010, 2011, 
2011a) with a symmetric scatter of data points around it, or it is 
caused by a boundary (as proposed by Nakar & Piran 2005), i.e. 
there is a considerable fraction of bursts with intermediateAarge 
Epcak and low luminosities (§5). We also verify with our code (§4) 
the recent claims on selection effects (Kocevski 2012) induced by 
the BATSE trigger threshold. 

The main advantage of our population synthesis code is that it 
relies on a small number of assumptions (described in §2). We focus 
on testing the Epeak — Liso correlation because for any simulated 
GRB with given z, Liso and Epeak, its 15-150 keV peak flux is 
easily compared with the Fum of the Swift complete sample. The 
adoption of the complete Swift sample is also relevant, because the 
Epcak — Liso correlation defined with this sample is independent 
from biases induced by the measurement of the redshifts (N12). 
Throughout the paper we assume a standard flat universe with h = 
^A = 0.7. 



2 SIMULATION SETUP 

We need to simulate a population of bursts with z and Liso as- 
suming a redshift density distribution and a luminosity function, 
respectively. To every simulated GRB we assign also a peak energy 
LJpeak and assume a typical spectral shape in order to compute its 
15-150 keV flux and compare it with Fuai- The main steps of our 
population synthesis code are: 

(1) We simulate a population of GRBs distributed in redshift z 
according to the GRB formation rate (GRBFR) ip{z) derived by Li 
2008 (which extended to higher redshifts the results of Hopkins & 
Beacom 2008): 

0.0157 + 0.118Z 



i>(z) 



(1) 



1+ (2/3.23)4-66 

where tp{z) is in units of Mq yr^^ Mpc"'^. This is the same ip{z) 
adopted by S 12 to derive the LF from the complete Swift sample. 
We simulate bursts with z ^ 10. 

(2) We assign peak luminosities to the simulated GRBs adopting 
the same luminosity function ^(Liso) used by S12 for studying the 
complete Swift sample: 



^iso ^ ^cut 



Liso > Lc 



(2) 



<;/)(Liso) oc ( --^^ 

Liso 
Lcut 

In S12 both the cases of a luminosity function evolving with red- 
shift as (1 + 2)*'" or a density evolution of the GRB population 
as (1 + 2) ° are considered. Therefore, we distinguish these two 
cases. Note that S12 assumed the -Epeak — Liso correlation in con- 
straining the LF free parameters (a, b, Lcut, ^l, 5d). Since this is 
in contrast with the null hypothesis we aim to test (i.e. that there is 
no Speak — Liso Correlation), we cannot use the parameter values 
of the LF reported in S12. For consistency with our null hypothe- 
sis, we have recomputed the LF of the complete Swift sample with 
the same analytic method adopted in S12 but without assuming the 
Speak — Liso coiTclation. The LF parameters (a, b, Lcut, (5l, 5d) 
adopted in our simulation (see Tab. [TJ are consistent, within the 
errors, with those reported in S12. 
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(3) We assign to each simulated burst a peak energy Speak- We 
consider that -Epcak is uncorrelated with Liso- We adopt a log- 
normal ^(-Bpeak) distribution centered at 337 keV (i.e. the average 
value for long GRBs with known z) with a dispersion of cr = 0.6 
dex. Assuming such a large dispersion of ^(Spcak) ensures that we 
simulate soft GRBs (i.e. X-ray flashes) with -Epoak of few tens of 
keV and hard bursts with -Epcak of tens of MeV. 

(4) We derive the peak flux of each simulated burst in the 15- 
150 keV energy band of Swift-BKi: 



F = 



N{E-a,p,E'X^,z)dE 



(3) 



where N{E; a, /3, iSpcJ^i ^) is the observer frame photon spectrum 
re-normalized through Liso. We assume for the ORB spectra the 
Band function (Band et al. 1993) with fixed low and high energy 
spectral index, a = —1.0 and /3 — —2.25 (as done in S12). 



From the simulated sample of GRBs we extract a sub-sample 
of bursts with (a) F ^Fii,„= 2.6 ph cm^^s"^ and (b) E°^|,, 
within 15 keV and 2 MeV. The GRBs used for the definition of the 
spectral-energy correlations have their Ep^ak measured. It is not 
mandatory that E°]^ly. is measured through Swift data, since most 
of the Swift bursts in the complete sample have been observed by 
the Konus Wind satellite or other instruments (such as the GBM on- 
board Fermi). Still we require that Bpeak is within the range of en- 
ergy covered by the present instruments. The sub-sample selected 
as above is the "complete sample of simulated GRBs". 

We then study the correlation between Epcak and Liso of the 
complete sample of simulated bursts and compare it with the ob- 
served correlation defined by the complete Swift sample. We will 
also perform the same simulation by modifying assumption (3). In 
particular we will assume that -Epcak and -Liso are coiTelated as 
found in real GRBs of the complete Swift sample (§5.1) or that the 
iSpcak — Liso coiTelation is caused by a boundary in the correspond- 
ing plane (§5.2). 



3 RESULTS 

By studying the complete Swift sample of GRBs, N12 found that 
the -Epcak — Eiso correlation defined by the 46 GRBs in this sample 
with known -Epcak and Liso has the following properties: 

• its rank correlation coefficient is p = 0.7 and its associated 
chance probability Pp = 10^^, i.e. it is significant at more than 
3a- 

• the fit of the correlation with a powerlaw log(-Epcak) = 
m log(Liso) + g gives m = 0.53 ± 0.06 (Icr) and g = -25.3 ±3.2 
(Icr) for the correlation slope and normalization, respectively; 

• the scatter of the 46 GRBs computed perpendicular around 
this best fit line has a dispersion of CTc — 0.29 dex; 

• the complete Swift sample has no outlier at more than 3cr, ex- 
cept for one burst at the limit of the 3a dispersion of the coiTelation. 

The goal of our simulations, for a given set of input assump- 
tions (i.e. tp{z), 0(Eiso), C(^peak), is (1) to producc a significant 
correlation in the Epcak — Eiso plane and (2) to obtain a coiTelation 
which is consistent with that of the real GRB sample of compari- 
son, i.e. the Swift complete sample. The simulation described in §2 
is repeated 300 times with the same initial assumptions. For each 
repeated simulation the complete sample of simulated GRBs is an- 
alyzed deriving the parameters of its Epcak — Eiso correlation. At 
the end of a cycle, i.e. after 300 repeated simulations, we derive: 
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Figure 1. Simulation of GRBs with density evolution (see Tab. [T]for the 
parameters). The simulated population of GRBs is shown by the dashed 
contour (1, 2, 3 and 4 a levels). The complete sample of simulated GRBs, 
i.e. those with a flux in the 15-150 keV energy range larger than the flux 
limit -Fiijn=2.6 ph cm~^ s~^ of the Swift complete sample, is shown by the 
(red) solid contours (1, 2 and 3 a levels). The open (blue) circles are the 46 
GRBs of the complete Swift sample analyzed in N12: the solid line is the 



best fit -E, 



peak 



Liso correlation to these bursts and the dot-dashed lines 



represent the 3(t scatter of the data points around the best fit line. 
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Figure 2. Simulation of GRBs with luminosity evolution (see Tab.[T]for the 
parameters). Same symbols as in Fig.[T] 



(i) the percentage Vp (Col. 6 in Tab.[Tll of simulations showing a 
£pcak~Eiso correlation of the complete sample of simulated GRBs 
significant at least at the 3a level. We consider a correlation signif- 
icant if the chance probability of its rank correlation coefficient is 
< 10'^ 

(ii) the percentage V (Col. 7 in Tab.[T]l of simulations significant 
at more than 3a with a slope m and normalization q consistent 
(within their \a errors) with the slope and normalization of the 
correlation defined by the Swift complete sample and with a scatter 
(Jc ^ 0.29, i.e. that of the observed correlation; 

(iii) the average percentage of simulated bursts which are out- 
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<J.(L) 


a 


-^cut 


erg s 1 


b 


<5 


-Pp 


V 


%Out. t 


%Out. i 


Density 

Luminosity 

Assume Ep<,ak - ^iso 

-E^pcak - Lisa boundary 


-1.37 
-L4 
-L4 
-1.4 


3.8x 

1051 
1051 


1052 


-2.37 
-2.13 
-2.13 
-2.13 


1.22 
2.67 
2.67 
2.67 


13% 
8.3% 
100% 
87% 


0.7% 
0.6% 
66% 
12% 


0.7% 
1.0% 
0.07% 

1.4% 


2.2% 
0.2% 
0.1% 


K12 (BATSE) 
K12 (Swift) 


-1.22 

-1.22 


1053 
1053 




-3.89 
-3.89 




0.5% 
0.7% 


0.0% 
0.0% 


2.6% 
0.4% 


0.7% 
1.3% 



Table 1. Vp is the percentage of simulations giving a significant coiTelation (i.e. with chance probability Sj 10~ ). V is the percentage of simulations giving 
a significant correlation with slope m and normalization q consistent, within their 1 a errors, with the correlation of the real Swift complete sample and with 
a scatter cfc ^ 0.29 (i.e. that of the observed coiTelation of the Swift complete sample). % Out J, (t) give the percentage of GRBs in the complete sample 
of simulated bursts which are outliers at more than 3(t of the i?peak ~ ^iso correlation defined by the complete Swift sample. The an'ows con'espond to the 
outliers below {\) and above (t) the boundary of the 3o" scatter (dot-dashed blue lines in all the figures). 



liers at more than 3(7 of the correlation defined by the Swift com- 
plete sample (Col. 8 and 9 in Tab.[TJ. 

For a certain set of assumptions {ij^^z) and (j)[Liso)) our term 
of comparison is the correlation observed in the Swift complete 
sample (N12). Therefore, the best match between the "simulated 
world" and the "reality" is when Vp and V are large and the per- 
centage of outliers at more than 3(j is consistent with 0.3% 

In addition, we also check that the complete sample of simu- 
lated GRBs has a redshift distribution and a flux distribution con- 
sistent with those of the real complete Swift sample. Since we adopt 
the LF derived from the complete sample of S12, we should find a 
complete sample of simulated bursts with ■(/'(z) and (j>{Liso) con- 
sistent with those of the real sample, and indeed this is the case. 

Tab.[T]lists the input assumptions (LF parameters) and the re- 
sults of our simulations which are also shown in Fig. [T] |2] |3] |4] 
|5] the contours are obtained by smoothing the distribution of data 
in the iSpcak — ^iso plane obtained by stacking all the 300 re- 
peated simulations. The red contours, for instance, are obtained by 
smoothing the distribution of 13800 data points, i.e. 46 simulated 
GRBs in 300 repeated simulations. 



NULL HYPOTHESIS: NO INTRINSIC £pcak 
CORRELATION 



.^is 



Under the null hypothesis that there is no intrinsic correlation be- 
tween E'pcak and Liso we find that ~92.7% (91.7%) of the 300 
repeated simulations do not produce a significant (at least at the 3(7 
level of confidence) Speak — ^iso correlation (Tab. [T) in the case 
of density (luminosity) evolution. Moreover, if we also require that 
the correlation obtained with the simulated bursts is consistent (in 
terms of its slope, normalization and scatter) with that observed in 
the Swift complete sample, the percentage of simulations satisfy- 
ing all our constraints drops to "P = 0.7% (V = 0.3% in the case 
of luminosity evolution). In other words, we can never reproduce 
the observed correlation through our simulations if we assume that 
Speak and Liso are uncorrelated. Therefore, we can reject the null 
hypothesis (i.e. that there is no correlation between Speak and Liso) 
at the 2.7(T and 3ct confidence level in the case of density and lumi- 
nosity evolution of the GRB population, respectively. These results 
are shown in Figs. [T]|2] where the solid contours represent the sub- 
sample of simulated bursts extracted with the same flux limit of the 
Swift complete sample (open blue circles in these figures) and the 
dashed contours represent the total population of simulated bursts. 
From these plots, it appears that the simulated bursts cannot repro- 



duce the distribution of the real data points (blue open circles) in 
the Speak — iiso plane. 

The solid (red) contours of Fig.[T|show a very weak correlation 
(in the upper-left part of the plane). Indeed, the flux-limit induces 
a weak correlation in the Speak — iiso plane by excluding part of 
the observables' space (i.e. the upper-left part of the Speak — Liso 
plane). However, a very weak correlation as that shown by the (red) 
solid lines in Fig. [T] has a very high significance if it is computed 
for all the 13800 data points resulting from the 300 repeated simu- 
lations. This is because the correlation significance would be dom- 
inated by the extremely large data sample. Instead, our goal is to 
compare samples of comparable sizes: the 46 GRBs of the swift 
complete sample versus the 46 GRBs of the simulated complete 
sample (the latter are produced in each of the 300 repeated simula- 
tions). For this reason, in analyzing our results, we have computed, 
for instance, the percentage of repeated simulations that generate a 
sample of 46 simulated GRBs which have a significant correlation 
in the Speak — Siso plane. As we discuss in §3, although an appar- 
ent very weak correlation is produced by the flux-limit in Fig.[T]in 
the uppper-left part of the plane, the overall distribution of the data 
points violates most of our constraints. 

Under the null hypothesis of no correlation between Speak 
and Liso , we have computed the average percentage of GRB s which 
are outliers at more than 3 a below (%Out J, in Tab.[T) the observed 
Speak — Liso Correlation defined by the Swift complete sample. Ac- 
cording to our simulations. Swift should have detected a percentage 
between 2.0-2.2% (see Tab. [TJ of GRBs with intermediate Speak 
and large Liso (e.g. lo52~54 gj.g ^-i) These bursts (1 event on av- 
erage) should be below the 3a limit of the Speak — Liso correlation 
(dot-dashed line in Fig.[T]|2]l defined by the complete Swift sample. 
At present, there are no outliers in the bottom/right triangle of the 
Speak — Liao plane, either in the complete Swift sample or more 
generally in the whole population of GRBs with measured z and 
known Speak and Liso (e.g. see N12). 

A similar argument applies to the simulated GRBs of the com- 
plete sample that have high Speak and low Liso- The 0.7-1.0% of 
simulated bursts with F ^Siim lie above the 3a upper boundary of 
the Speak — Liso Correlation of the real burst sample. While there 
are no outliers in the complete Swift sample in this part of the plane, 
it has been discussed in the literature if a fraction of bursts with in- 
termediate/large Speak and low Liso could already be present in the 
populations of GRBs detected by different satellites (e.g Nakar & 
Piran 2005; Band & Preece 2006). On the other hand, Nava et al. 
(2011) have shown that there are no outliers of the Speak — Liso 
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Figure 3. Results of the simulation of a GRB population with (f){L[ao) and 
^(z) derived by Butler et al. 2011 (used by K12) and adopting the BATSE 
trigger threshold (solid contours) or the Swift i<iini(dot-dashed contours). 
The null hypothesis is that there is no intrinsic iJpeak — ^iso correlation. 
Symbols as Fig.[T] 



correlation at more than 3 it in the BATSE and Fermi population of 
GRBs. 

Note also that a considerable fraction (~20%) of the simulated 
GRBs with F ^Fu^ have a peak energy SpeJ],., in the observer rest 
frame, larger than the upper Swift energy threshold of 350 keV but 
still lower than 2 MeV, which roughly corresponds to the limit of 
current instruments with the largest energy range that can measure 
^pclk (s-g- Fermi and Komis). 

Kocevsky 2012 (K12, hereafter) performed a population study 
finding that the Speak — ^iso correlation is induced by a combina- 
tion of the redshift/luminosity function of GRBs with the detection 
limit of a given instrument (BATSE in his study). K12 attributes 
the existence of the i?peak — -Eiso correlation to the combination of 
the Malmquist bias (i.e. the detection of the most luminous GRBs 
at high z preferentially), the flux limit of the detector, and its lim- 
ited band-pass. According to this interpretation the most luminous 
GRBs at high z with intermediate -Epoak goes undetected because 
their i5peJij=_Epeak/(l + z) falls below the low energy threshold of 
the detector (i.e. 15 keV in the case of Swift-UPS). This argument 
is intuitively plausible: a detector sensitive in a given energy range 
should hardly detect bursts with Bp^ak outside this range. However, 
what matters for triggering a burst is its peak flux (integrated in 
the detector energy range) which depends on the spectral shape of 
the burst (e.g. -Bpcak^ ct> P) and the normalization of the spectrum. 
If this flux is larger than the detector flux limit, the burst will be 



detected no matter its E\ 



peak 



Our simulations in the case of lumi- 



nosity evolution, for instance, predict that ~13% of the simulated 
GRBs with F ^Fu,^ have £°J^|k<15 keV (i.e. the lower thresh- 
old of the BAT energy band-pass). In other words, these bursts are 
bright enough to be detected by Swift, despite the fact that their 
^peak is outside the BAT energy range. Therefore, the absence of 
GRBs in the lower/right plane of the Bpcak — iiso correlation can- 
not be entirely attributed to this effect. 

Our population synthesis code is simpler than that of K12 be- 
cause we test the Speak — Siso correlation through the complete 
Swift sample. In our case, by simulating z, Liso and Speak we 
can immediately compute the peak flux of a simulated GRB and 



compare it with the flux limit. The simulation of K12 concerns the 
Speak — Siao Correlation, which involves the time integrated spec- 
tral properties of GRBs, and requires several assumptions (e.g. on 
the profile structure and time evolution of the spectrum during a 
burst) in order to simulate a GRB light curve and verify if it can 
be detected by BATSE (in turn, this also requires K12 to model a 
typical background and the BATSE detector response matrix). The 
advantage of our simulation is that it uses the Sum of the complete 
sample of S 12, which is a simple sharp cut in the 15-150 keV peak 
flux of the simulated GRB population. 

K12 adopts a slightly different LP and performs the simulation 
for BATSE. We assumed his luminosity function and redshift den- 
sity distribution <j}(Liso) and ^(z) (originally derived from Butler et 
al. 201 1). We have implemented in our code the BATSE detection 
algorithm as described in Band (2003). The results of this simula- 
tion are shown in Fig.|3]and reported in Tab. [T|(labelled K12). Un- 
der the null hypothesis of no correlation between Speak and Liso, 
the majority (99.5%) of the repeated simulations show no signifi- 
cant Speak — Liso Correlation in the sample of bursts that should be 
detected by BATSE. The percentage of outliers that we find is 0.7% 
below the correlation and 2.3% above the correlation (see Fig. [3] 
and Tab.[TJ. We compare the results of the simulation assuming the 
BATSE trigger threshold with the complete Swift sample instead 
of using the larger (but highly incomplete) sample of all the bursts 
with measured z and known spectral parameters. This is justified by 
the findings of N12 that show that the Speak — Liso correlation de- 
fined by the complete Swift sample (46 events with Speak and Liso 
known) is quite similar to the correlation defined with the larger 
sample of 136 GRBs with measured z and known spectral parame- 
ters (Speak, Liso)- We have also tested the LP adopted by K12 but 
with the flux limit of Swift. We show in Pig.[3](see also Tab.lTJ the 
contours coiTesponding to the complete sample of simulated GRBs 
(dot-dashed contours). Also in this case we are not able to obtain 
a significant (at more than 3a") Speak — Liso correlation in 96% of 
the repeated simulations under the null hypothesis that there is no 
Speak — Liso Correlation. 



5 S, 



peak 



Liso CORRELATION OR BOUNDARY? 



We have shown in §4 that if there is no correlation between Speak 
and Liso the apparent correlation between these two observables 
cannot be produced by the flux limit cut of the complete Swift sam- 
ple or of the BATSE instrument. We still have two possibilities: 
(i) to assume the observed correlation and use it in simulating the 
GRB population; (ii) assume that the Speak — Liso correlation is 
produced by a boundary in the corresponding plane. 

Por these reasons, it is worthwhile to test the possibility that 
there is a population of GRBs with a uniform distribution of Speak 
above the Speak — Liso correlation. In this case their absence in the 
observed Speak — Liso plane could be induced by the trigger thresh- 
old and/or by a bias related to the detection of the bursts with the 
smallest jet opening angles. Here we study the former possibility 
(i.e. absence due to the detector threshold). 



5.1 Real correlation 

We repeat the Monte Carlo simulation to obtain a population of 
GRBs assuming that Speak and Liso are intrinsically correlated. We 
assume that the correlation between Speak and Liso has the same 
properties as the correlation that we observe in the complete sample 
of Swift bursts (point (3) of §2). In particular we refer to the slope. 
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Figure 4. Simulation of GRBs with luminosity evolution (see Tab. [T] for 
the parameters) assuming the real Speak ~ ^iso correlation observed in 
the complete Swift sample of GRBs. Same symbols as in Fig.[T] 



Figure 5. Simulation of GRBs with luminosity evolution (see Tab.[T]for the 
parameters) assuming the real -Epeak ~ ^iso correlation as a boundary: the 
distribution probability of -Epcak is from Eq.|4] Same symbols as in Fig.[T] 



normalization and scatter found in N12. Then, from our simulated 
population of GRBs, we select those that satisfy the requirements 
on the peak flux and on -Bpeak and we study the Speak — iiso 
correlation in this sub-sample. In this sub-sample we expect, of 
course, to find a strong correlation, since we have assumed it in 
simulating the population of GRBs. 

The main reason for performing this test is for self- 
consistency and also to verify wether the cuts introduced (flux limit 
and constraints on -Bpelk) in defining the complete sample of simu- 
lated GRBs do affect the characteristics of the correlation and lead 
us to find a correlation which is different from the correlation that 
we have assumed for the simulation. For example, if we simulate 
a sample satisfying a given correlation, we can expect that the in- 
troduction of a flux limit could reduce the scatter (and maybe also 
modify the slope) of the correlation that we find. 

The results are shown in Fig.|4]and reported in Tab. [T] We find 
a good agreement between the complete Swift sample of simulated 
bursts (solid contours in Fig. |4j and the complete Swift sample of 
real bursts (open points in Fig.|4j. All the repeated simulations (we 
have assumed the LF with luminosity evolution but the results do 
not depend on this choice) produce a significant correlation in the 
Speak — Lisa plane (see Tab.[TJ and in 66% of the simulations such 
a significant correlation is consistent in slope, normalization and 
scatter with that defined by the Swift complete sample. 



5.2 Boundary 

We can also test the possibility that the Speak — Siso correlation 
is caused by a boundary. To this aim, we peiform the simulation 
by assuming that bursts follow the boundary represented by the 
log(Speak) ~ mlog(I/iso) + Q (with q and m being the normal- 
ization and slope of the Speak — Siso correlation found by N12 for 
the complete Swift sample). For each simulated value of Liso, we 
assign a peak energy Speak according to a probability distribution: 



<^(Spe 



Ae~ 

A; 



-C^peak ^ -L/peak,c 
^peak ^ ^peak,c ('+) 



where log(Spcak,c) ~ n^log(Liso) + ?• This distribution intro- 
duces a sharp decrease of the simulated bursts below the observed 
Speak ^ Lisa Correlation, which then forms a boundary in the 
Speak — Liso plane with a uniform distribution of GRBs above 
it. We report the results in Tab. [T]and show them in Fig. |5] Under 
this hypothesis, we find that a considerably large fraction (87%) of 
the repeated simulations produces a strong Speak — Siso correlation 
but only 12% of the correlations seem consistent with the correla- 
tion of the Swift complete sample. 1.4% of outliers at more than 
3a above the correlation are found. By contrast, in the samples of 
bursts without measured redshifts of BATSE and Fermi, Nava et al. 
(201 1) have shown that there are no outliers at more than 3 a of the 
Speak — Siso correlation. 



6 CONCLUSIONS 

We performed Monte Carlo simulations of GRBs with a given red- 
shift distribution and luminosity function (as described in Eq. [T] 
and Eq. |2] respectively). We have made the null hypothesis that 
there is no intrinsic correlation between the luminosity Liso and 
the peak energy Speak of the simulated bursts. We have consid- 
ered the GRBs that have a peak flux Js Siim= 2.6 ph cm^'^ s~^ 
which is the same flux limit of the complete Swift sample studied 
in S12 and N12 and we also required that E°]^l^^ of the detected 
bursts can be measured by current instruments, i.e. that it lies in the 
15 keV-2 MeV energy range. These are the simulated bursts that 
would be detected by Swift. The use of a flux-limited sample of 
Swift bursts as that defined in S12 has the advantage that it avoids 
several complexities related to the Swift trigger method (e.g. depen- 
dence from the time evolution of the signal and significance with 
respect to a time-variable background or dependence from off-axis 
detector response). In particular the flux-limit of the S 12 sample is 
sufficiently high that the selected sample should be free from these 
trigger-related issues. 

If we make the hypothesis that there is no correlation between 
Speak and Liso, only in 7.3% of the repeated simulations (e.g. for 
the case of a GRB population evolving in density with redshift, as 
found in S12) we find a statistically robust (i.e. chance probability 
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of the rank correlation coefficient ^ 10~^) Spoak ^ ^iso correla- 
tion. If we also require that our simulations produce a correlation 
similar (in slope, normalization and scatter) to that observed among 
real GRBs of the complete Swift sample, the percentage reduces to 
0.7%. Therefore, we reject the null hypothesis that there is not an 
intrinsic Speak — iiso correlation at the 2.1a level of confidence 
(3.0(7 for the case of luminosity evolution). These results suggest 
that a correlation between Liso and Speak should exist. Since for 
Swift we have considered a bright cut on the peak flux (to match 
that adopted in the definition of the complete sample of S12), our 
results are obtained for the most conservative case. Similar results, 
i.e. the impossibility to produce a strong -Epeak ^ iiao correlation 
as due to the trigger threshold, are also obtained for BATSE. 

An alternative possibility is that there is a boundary in the 
Speak — -^iso plane (Nakar & Piran 2005). If we assume that this 
boundary is coincident with the observed Speak — iiso correlation, 
the scatter of the data points around it is not symmetric but there is 
a larger fraction of bursts with intermediate Speak and low lumi- 
nosities. The trigger threshold here could play a role in hiding the 
events with extremely low luminosity. Our simulations, assuming 
the existence of the Speak — iiso correlation as a boundary, suggest 
that while a considerable percentage (87%) of simulations produce 
a significant correlation in the Speak — Lisa plane, only 12% of 
the simulations also produce a correlation with slope, normaliza- 
tion and scatter consistent with that of the real GRB sample. This 
result excludes at the 1. 6a confidence level the boundary case, as 
we have modeled it here. In particular, we assumed that the dis- 
tribution of bursts above the boundary is uniform. It could still be 
possible, however, that a boundary with an asymmetric scatter be- 
low and above it exists due to some other property of GRBs (e.g. 
their jet opening angle). This possibility requires more elaborated 
simulations that are outside the scope of this paper In conclusion, 
our simulations show that there is a correlation between Speak and 
Liso and this caimot be due to a selection bias caused by a flux- 
limited sample. 
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